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Abstract. The universe cools down monotonically following its expansion. This generates a se- 
quence of phase transitions. If a second order phase transition happens during the radiation domi- 
nated era with a charged order parameter, spinodal unstabilities generate large numbers of charged 
particles. These particles hence produce magnetic fields. We use out of equilibrium field theory 
methods to study the dynamics in a mean field or large N setup. The dynamics after the transition 
features two distinct stages: a spinodal regime dominated by linear long-wavelength instabilities, 
and a scaling stage in which the non-linearities and backreaction of the scalar fields are dominant. 
This second stage describes the growth of horizon sized domains. We implement a formulation 
based on the non-equilibrium Schwinger-Dyson equations to obtain the spectrum of magnetic fields 
that includes the dissipative effects of the plasma. We find that large scale magnetogenesis is ef- 
ficient during the scaling regime. Charged scalar field fluctuations with wavelengths of the order 
of the Hubble radius induce large scale magnetogenesis via loop effects. The leading processes are: 
pair production, pair annihilation and low energy bremsstrahlung, these processes while forbidden in 
equilibrium are allowed strongly out of equilibrium. The ratio between the energy density on scales 
larger than L and that in the background radiation r(L, T) — Pb(L, T) / Pcmb{T) is r(L 1 T) ~ 10~ 34 
at the Electroweak scale and r(L, T) ~ 10~ 14 at the QCD scale for L ~ 1 Mpc. The resulting spec- 
trum is insensitive to the magnetic diffusion length and equipartition between electric and magnetic 
fields does not hold. We conjecture that a similar mechanism could be operative after the QCD 
chiral phase transition. 
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EARLY COSMOLOGY AND FUNDAMENTAL PHYSICS 

The history of the universe is determined by its expansion and consequent cooling. 
During most of its early history the Universe is homogeneous and isotropic to an 
excellent approximation and is therefore described by the spatially flat Friedmann- 
Robertson-Walker (FRW) geometry 

ds 2 = dt 2 -a 2 (t) dx 2 (1) 
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where the scale factor a(t) grows with t. Physical lengths are proportional to a(t) and 
the temperature decreases as T(t) ~ ^Ly. This monotonous decrease of the temperature 
generates a sequence of cosmological phase transitions with the ensuing breaking 
of internal symmetries. The symmetry of the Universe reduces through each phase 
transition. 

The main ingredients to describe the early Universe are: 

• General Relativity: Einstein's Theory of Gravity 
The matter distribution determines the geometry of the spacetime through the 
Einstein equations. For the geometry eq.©, the Einstein equations reduce to one 
scalar equation, the Einstein-Friedman equation 



da 



a(t) dt 



^Gp(t), (2) 



where G stands for Newton's gravitational constant and p(t) for the energy density. 

Quantum Field Theory and String Theory to describe Matter 

Since the energy scale in the early universe is so high (well beyond the rest mass 

of particles), a quantum field theoretical description for matter is unavoidable. 

Only such context permits a correct description of particle production and particle 

decays. 



Electromagnetic, weak and strong interactions are well described by the so-called 
standard model. That is, quantum chromodynamics (QCD) combined with the elec- 
troweak theory (electromagnetic and weak interactions). This a non-abelian gauge the- 
ory associated to the symmetry group SU (3) £g> SU (2) ® U ( 1 ) . The SU (3) corresponding 
to the color group of QCD while SU (2) <g> U(l) describes the electroweak sector. To this 
scheme, one adds presently neutrino masses (through the see-saw mechanism) to explain 
neutrino oscillations. 

The energy scale in QCD is about ~ lOOMeV ~ 10 12 K corresponding to the chiral 
symmetry breaking and determined by the pion mass, while the energy scale for the 
electroweak is the Fermi scale ~ lOOGeV ~ 10 15 K, which is determined by the mass of 
the vector bosons. 

The standard model has been verified experimentally with spectacular precision. 
However, it is an incomplete quantum field theory and it is the major challenge of our 
times to understand its extension. It seems obvious that extensions of the Standard model 
will be symmetric under a group containing SU (3) <E> SU (2) ® U(l) as a subgroup. Pro- 
posals for a Grand Unified Theory (GUT) include 50(10), SU(6) and £g as symmetry 
group. 

The grand unification idea consists in that at some energy scale all three couplings 
(electromagnetic, weak and strong) should become of the same strength. The running 
of the couplings with the energy (or the length) is governed by the renormalization 
group. For the standard model of electromagnetic, weak and strong interactions, the 
renormalization group yields that the three couplings get unified approximately at ~ 
10 16 GeV yj]. A better convergence is obtained in supersymmetric extensions of the 
standard model|[l|]. 



Grand unified models may posses magnetic monopoles or not according to their 
symmetry group and to the symmetry of the ground state. Notice that no experimental 
evidence for magnetic monopoles has been found so far. 

Quite generally, the internal symmetry increases with energy. This is true in general, 
in statistical mechanics, condensed matter as well as in cosmology. For example, a 
ferromagnet at temperatures higher than the Curie point is in the symmetric phase 
with zero magnetization. Below the Curie point, the internal symmetry is spontaneously 
broken by a non-zero spontaneous magnetization. 

Current models purport that the universe started with maximal symmetry before 
inflation and this symmetry reduces gradually while the universe expands and cools. 
The symmetry breaking transitions includes both the internal symmetry groups (as the 
GUT's symmetry group that eventually reduces to the SU(3) (g>Sl/(2) ® E/(l) group) 
as well as the translational and rotational symmetries which are broken by the density 
fluctuations amplified by gravitational instabilities leading to structure formation. 

It should be noticed, however, that no direct manifestation of supersymmetry is known 
so far. An indication emerges by studying the energy running of the (electromagnetic, 
weak and strong) in the standard model and in its minimal supersymmetric extension 
(MSSM). All three couplings meet at E ~ 2 x 10 16 GeV in the MSSM. The coupling 
unification becomes quite loose in the Standard Model. This is why the renormalization 
group running of the couplings in the MSSM supports the idea that supersymmetry 
would be a necessary ingredient of a GUT. 

Neutrino oscillations and neutrino masses are currently explained in the see-saw 
mechanism as follows!!, 



M 

where Mf erm j ~ 250 GeV is the Fermi mass scale, M ^> Mp erm i is a large energy scale and 
Am v is the difference between the neutrino masses for different flavors. The observed 
values for Am v ~ 0.009 — 0.05 eV naturally call for a mass scale M ~ 10 15 ~ 16 GeV close 
to the GUT scaled. 

The third evidence for an energy scale about 10 16 GeV comes from inflation: data on 
the cosmic microwave background (CMB) anisotropics indicate that the inflation scale, 
the grand unification scale and the supersymmetry breaking scale actually coincide |0]. 



PRIMORDIAL SEEDS FOR THE MAGNETIC FIELDS IN THE 

UNIVERSE 

A variety of astrophysical observations including Zeeman splitting, synchrotron emis- 
sion, Faraday rotation measurements (RM) combined with pulsar dispersion measure- 
ments (DM) and polarization measurements suggest the presence of large scale magnetic 
fieldso,0]. The strength of typical galactic magnetic fields is measured to be ~ \i G|0,0] 
and they are correlated on very large scales up to galactic or even larger reaching to 
scales of cluster of galaxies ~ 1 Mpc[5]. The origin of these large scale magnetic fields 
is still a subject of much discussion and controversy. It is currently agreed that a vari- 
ety of dynamo mechanisms are efficient in amplifying seed magnetic fields with typical 



growth rates T ~ Gyr over time scales ~ 10—12 Gyr (for a thorough discussion of 
the mechanisms and models seeo]). The ratio of the energy density of the seed mag- 
netic fields on scales larger than L (today) to that in the cosmic background radiation, 
r(L) = Pb(L) I pcmb must be r(L ~ IMpc) > 10~ 34 for a dynamo mechanism to amplify 
it to the observed value, or r(L ~ IMpc) > 10~ 8 for the seed to be amplified solely by 
the gravitational collapse of a protogalaxy[5]. 

There are also different proposals to explain the origin of the initial seed. Astrophysi- 
cal batteries rely on gradients of the charge density concentration and pressure and their 
efficiency in producing seeds of the necessary amplitude is still very much discussed^. 
Primordial magnetic fields that could be the seeds for dynamo amplification can be gen- 
erated at different stages in the history of the early Universe, in particular during in- 
flation, preheating and or phase transitions! 5]. Primordial (hyper) magnetic fields may 
have important consequences in electroweak baryogenesis[7], Big Bang nucleosynthesis 
(see^), the polarization of the CMBli via the same physical processes as Faraday rota- 
tion, and structure formation^ Ejj], thus sparking an intense program to study the origin 
and consequences of the generation of magnetic fields in the early Universe! 101- fT^l . 

A reliable estimate of the amplitude and correlations of seed magnetic fields must 
include the dissipative properties of the plasma, in particular the conductivity ifTli foil . 
In ref . II 1 9fl we have introduced a formulation that allows to compute the generation of 
magnetic fields from processes strongly out of equilibrium. This formulation, which 
is based on the exact set of Schwinger-Dyson equations for the transverse photon 
propagator is manifestly gauge invariant and is general for any matter fields and any 
cosmological background (conformally related to Minkowski space-time). In the case 
in which strongly out of equilibrium effects arise from long-wavelength fluctuations, 
such as during phase transitions, this formulation allows to separate the contribution of 
the hard degrees of freedom which are in local thermodynamic equilibrium from that 
of the soft degrees of freedom that fall out of LTE (local thermal equilibrium) during 
the phase transition and whose dynamics is strongly out of equilibrium. This separation 
of degrees of freedom leads to a consistent incorporation of the dissipative effects via 
the conductivity (for details see[19]). In that reference a study of magnetogenesis in 
Minkowski space-time during a supercooled phase transitions was presented and the 
results highlighted the main aspects of the generation of magnetic and electric fields in 
these situations. 

We study the generation of large scale (hyper) magnetic fields by a cosmological 
phase transition during a radiation dominated era by implementing the formulation 
introduced in ref. lfl9ll . The setting is a theory of N charged scalar fields coupled to an 
abelian gauge field (hypercharge). We consider the situation when this theory undergoes 
a phase transition after the reheating stage and before either the Electroweak or the 
QCD phase transition, since we expect that these transitions will lead to new physical 
phenomena. The non-perturbative dynamics out of equilibrium is studied in the limit of 
a large number N of (hyper) charged fields and to leading order in the gauge coupling. 
The non-equilibrium dynamics of the charged scalar sector features two distinct stages. 
The first one describes the early and intermediate time regime and is dominated by 
the spinodal instabilities which are the hallmark of the process of phase separation 
and domain formation and growth. This stage describes the dynamics between the 
time at which the phase transition takes place and that at which non-linearities become 



important via the backreaction. The second stage corresponds to a scaling regime which 
describes the slower non-equilibrium evolution of Goldstone bosons and the process 
of phase ordering! 18] and growth of horizon-sized domains. This scaling regime is 
akin to the solution found in the classical evolution of scalar field models with broken 
continuous symmetries after the phase transition that form the basis for models of 
structure formation based on topological defectsl20. l2lll . 

The solution of the scalar field dynamics HI 8fl is the input in the expression for the 
spectrum of the magnetic field obtained in lll9ll to obtain the amplitude of the primordial 
seed generated during both stages. 

We find that scaling stage is the most important for the generation of large scale mag- 
netic fields. Large scale magnetic fields are generated via loop effects from the dynamics 
of modes that are at the scale of the horizon or smaller. The leading order processes that 
result in the generation of large scale magnetic fields are: i) pair production, ii) pair 
annihilation and iii) low energy bremsstrahlung. These processes would be forbidden 
in equilibrium by energy momentum conservation, but they are allowed strongly out of 
equilibrium because of the rapid time evolution of the cosmological background and the 
fast dynamics of the scalar field fluctuations. 

The resulting spectrum is rather insensitive to the diffusion length scale which is 
much smaller than the horizon during the radiation dominated era. The ratio of the 
magnetic energy density on scales larger than L (today) to the energy density in the 
background radiation r(L,r\) = ps(L, r\)/ p cmb {r\) is summarized in a compact formula 
req.(E0b1. For L ~ 1 Mpc (today) we find r(L, 77) ~ 10~ 34 at the Electroweak scale and 
r(L,7]) ~ 10~ 14 at the QCD scale, suggesting the possibility that these primordial seeds 
could be amplified by dynamo mechanisms to the values of the magnetic fields consistent 
with the observed ones on these scales. 

THE PHYSICAL PICTURE 

The extreme energy scale and energy density during the inflationary and radiation 
dominated eras call for a quantum field theoretic treatment of the matter and radiation 
while the geometry is described by the classical metric eq.©. The fast expansion 
of the universe can lead to out of thermal equilibrium situations, which require the 
implementation of out of equilibrium methods. In addition, nonperturbative methods 
are needed since the energy density is proportional to the inverse of the coupling. We 
developed nonperturbative field theory methods that successfully treats the inflationary 
era in various relevant scenarios and allowed to compute the primordial perturbations 
and the CMB fluctuations as well as to make contact with the customary slow-roll 
classical treatment I17I1 . 

We consider here a scalar field carrying an abelian charge and coupled to the electro- 
magnetic field in the radiation dominated era. During the second order phase transition 
the concavity of the potential becomes negative at the origin and strong spinodal fluctu- 
ations are generated. These fluctuations in turn generate a magnetic field with a typical 
wavelength corresponding to the wavelength of the spinodally unstable modes. 

This is the main premise of our work lll9ll : the instabilities which are the hallmark of 
a non-equilibrium symmetry breaking phase transition lead to strong fluctuations of the 



charged scalar fields which in turn, lead to the generation of magnetic fields through the 
non-equilibrium evolution. 

The main ingredients developed in ref. [19] to compute the generation of magnetic 
fields through this non-equilibrium process were: 

• A consistent framework to compute the spectrum of generated magnetic field, 
namely (B(k,t) ■ B(—k,t)) /V with B{k,t) the spatial Fourier transform of the 
Heisenberg magnetic field operator and V the (comoving) volume of the system. 

• Plasma effects were included to assess the generation and eventual decay of the 
magnetic fields. For a large conductivity in the medium, the magnetic field diffuses 
but also its generation is hindered. This point is of particular importance within the 
cosmological setting 1I5I H3L [Till. 

• A major challenge of any mechanism of large scale magnetogenesis is to generate 
the seed magnetic fields from microscopic, causal processes. An important aspect 
of the results presented here is that this generation mechanism is mediated by loop 
effects and correspond to processes that are forbidden in equilibrium but allowed 
strongly out of equilibrium. 



FIELD THEORETICAL MODEL FOR MAGNETIC FIELDS IN 
FRIEDMANN-ROBERTSON- WALKER COSMOLOGY 

We will not attempt to study a particular gauge theory phenomenologically motivated 
by some GUT scenario, but will focus our study on a generic scalar field model in 
which the scalar fields carry an abelian charge. The simplest realization of such model 
is scalar electrodynamics with N charged scalar fields <j> r , r = 1, . . . , N and one neutral 
scalar field \y whose expectation value is the order parameter associated with the phase 
transition. The neutral field is not coupled to the gauge field and its acquiring an 
expectation value does not break the U(l) gauged symmetry. This guarantees that the 
abelian gauge symmetry identified with either hypercharger or electromagnetism is not 
spontaneously broken to describe the correct low energy sector with unbroken U{V)em- 
We will take the neutral and the N complex (charged) fields to form a scalar multiplet 
under an 0(2N + 1) isospin symmetry. The electromagnetic coupling explicitly breaks 
the 0(2N + 1) symmetry down to SU(N) x 2/(1) . In the absence of electromagnetic 
coupling as the neutral field acquires an expectation value the isospin symmetry is 
spontaneously broken to 0(2N). Since by construction only the neutral field acquires 
a non vanishing expectation value under the isospin symmetry breaking the photon 
remains massless (it will obtain a Debye screening mass from medium effects. 

The action that describes this theory in a general cosmological background and using 
conformally rescaled fields is given by 

S = Jdri d 3 x ^d^d^ + D^*D^^-M 2 (ri) (*!L + &*&J 



AN V 2 ) 4 MV 



(3) 



with 



M z (r]) = ~^ 2 C 2 (ri)- 



c"(v) 



D, 



ieA, 



d^A v - d v A^ ,(4) 



and the primes refer to derivatives with respect to conformal time. Obviously the confor- 
mal rescaling of the metric and fields turned the action into that of a charged scalar field 
interacting with a gauge field in flat Minkowski space-time, but the scalar field acquires 
a time dependent mass term 2 . In particular, in the absence of electromagnetic coupling, 
the equations of motion for the gauge field A^ are those of a free field in flat space 
time. This is the statement that gauge fields are conformally coupled to gravity and no 
generation of electromagnetic fields can occur from gravitational expansion alone with- 
out coupling to other fields or breaking the conformal invariance of the gauge sector. 
The generation of electromagnetic fields must arise from a coupling to other fields that 
are not conformally coupled to gravity, or by adding extra terms in the Lagrangian that 
would break the conformal invariance of the gauge fields II 1 lh . 

The dynamics is determined by the Heisenberg equations of motion of the neutral field 
and the charged fields 4> HITL llsL llQII . We will consider that at the onset of the radiation 
dominated era, the system is in the symmetric high temperature phase in local thermal 
equilibrium with a vanishing expectation value for the scalar fields. In the absence of 
explicit symmetry breaking perturbations the expectation value of the scalar field will 
remain zero throughout the evolution, thus <p = 0. 

It is convenient to introduce the mode expansion of the charged fields 



d 3 k 



flr(*) Aft) + 



—ik-x 



r=l,...,N 



(5) 

In leading order in the large N limit, the Heisenberg equations of motion for the charged 
fields translate into the following equations of motion for the mode functions for r\ > 



^ + k 2 -M 2 (r 1 ) + ^cp 2 (r 1 ) + ^(&<i>) 



A(*7) = o. 



With our choice of the initial state we find the backreaction term to be given by| 1 



v u 2N \ i 4 



d 5 k 



\f k {r))\ 2 [\+2n k 



(6) 



(7) 



This expectation value features quadratic and logarithmic UV divergences which are 
absorbed in the mass and coupling renormalization ll 1 811 . 

After renormalization and in terms of dimensionless quantities, the non-equilibrium 
dynamics of the charged scalar fields is determined by |I17LI18LI19I1 . 



dr\ : 



/«(»?) =o ; 



(8) 



Here we neglect the effect of the conformal anomalv lll5ll 



with the effective time dependent mass given by 



T 2 = 



C 2 (t])M 2 

24 ju 2 
X+3 



C 2 (t?) T 2 



- 1 



,2 • 



(9) 



(10) 



We see that Jt 2 (r\) > for early times when C(tj) < T/T c . For later times C(tj) > T/T c 
and ^# 2 (r/) becomes negative triggering the phase transition at a time r/ c when C(r\ c ) = 
T/T c since we can neglect the nonlinear contribution AE(tj) (recall that A 1). 

When M 2 {r\) is negative, eq.© tells us that the modes / 9 (t]) grow exponentially 
aslO 



A r?>l 



fq(l) = a q e? 



>(£ I?) 



CM 



— 9 —i 



l 



fi 2 r\ 2 



(11) 



where jx = tJJTHr and is the Hubble constant at the reheating time, = r\ R 1 . 

This growth continues till the nonlinear term in eq.® E(tj) = (<J>^<I>) cannot be ne- 
glected anymore and stops the unstabilities. An scale invariant stage follows 111 811 . In this 
scaling stage the tree level mass term —C 2 {r\) jU 2 is compensated by the backreaction of 
the quantum fluctuations as follows 111 8(1 . 



AE(7 7 )- At 2 C 2 (7 ] ) r? = 



15 



(12) 



The mode functions can be written during the scaling stage in terms of Bessel functions 



m)=A k ^/ 2j ^ + B. 



x 2 N 2 {x) 



3/2 



(13) 



where x = k r\ is the scaling variable. The correlation length of the scalar field is of the 
order of the Hubble radius in this stage. 



Gauge field dynamics 

The electric conductivity is very large in the high temperature plasma formed after 
the second order phase transition and dominates the dynamics of the gauge field. The 
conductivity is obtained from the imaginary part of the photon polarization and it is 
dominated by charged particles of momenta p ~ T in the loop with exchange of photons 
of momenta eT < k <^ T. The effect of Debye (electric) and dynamical (magnetic) 
screening via Landau damping is crucial leading to the expression ll22ll . 

„(„)=* , « = ™ (H) 
C(tj) ahi^ 

where = 15.698 Such large conductivity leads to dissipative processes which 

severely hinders magnetogenesis (see lfl9ll ) and also introduces the diffusion length scale 
which could limit the correlation of the magnetic fields that are generated. 



As a consequence, we found for the photon causal correlator for k <^ Or 



® R {r),r)';k) = e{T)-T)') 



e -7r R (.n-i}') 
Or 



(15) 



MAGNETIC FIELD SPECTRUM 



The magnetic energy at wavenumber k is given by the symmetric equal time limit 



S B (rj,k) = \ Mm / d 3 x< {^(rj,!),^' ,0)} > p e 

l 77'— »T7 J 



ik-x 



(16) 



where { , } denotes the anti commutator, B(r),x) above is a Heisenberg operator and 
the expectation value is in the initial density matrix. 

The physical magnetic energy density stored on comoving length scales larger than a 
given L is given by 



1 ~ 

ApB {L, 7] ) = Jq L k 2 S B (ri,k) dk . 



(17) 



Apg(L, rj) stands for the contribution from the non-equilibrium generation (subtracting 
the local thermodynamic equilibrium contribution), a quantity of cosmological relevance 
to assess the relative strength of the generated magnetic field is given by the ratio of the 
power on scales larger than L to the energy density in the radiation background 



Ap g (M) 

Pr(n) 



Py 



til 
15 



(18) 



The explicit field theoretic evaluation of 5g(r],^) is given in ref.[19]. The leading 
contribution to the power spectrum generated by non-equilibrium fluctuations results 
expressed in terms of the mode functions fk(f]) as follows, 



S B {ri,k) = {l+2n ) 



aNk 2 



7t On 



e °r 



q A dq d(cos6) (1 -cos 2 0) 



J R 6 ° R f^ftf+kiMdm 



(19) 



where 6 is the angle between the vectors q and k and we have used the result that both the 
spinodal stage as well as the scaling stage is dominated by the long-wavelength modes 
that acquire non-perturbatively large amplitudes 111 911 . 

We then compute r(L, 17) using the scaling solution eq. (fT3t for the mode functions. 
The result can be recasted for L 3> Tj as lll9ll . 



r(T],L) = 3.665 x 10 4 



a 
Ob 



1 3 



1 



[LTr]' 



(20) 



where 



t r 



1 



M, 



Pi 



15.698... 



1^ 



C(V) VG 
Several important features of the above result are noteworthy: 

• i: A typical nonperturbative behaviour p- arising from the phase transition. 

• ii: r(r],L) grows with the symmetry breaking scale \x as /i 4 since the higher is \l 
the longer is the scaling stage. 

• iii: The presence of the huge suppression factor {LTr]^ 5 . 

For the symmetry breaking scale /i ~ 10 13 Gev and X ~ a ~ 10~ 2 , corresponding 
to a critical temperature of order of a GUT scale T c ~ 10 15 GeV and assuming that 
the scaling regime lasts until the electroweak phase transition scale, i.e. r\ is such that 
T{r\) = T EW ~ 10 2 GeV. Then the factor 



1" 



EW 



Tew 



10 



100 



compensates for the factor [LTr] 5 . Taking N and g* of the order of 10 we obtain 



r[T(rj),L}~ 10- 



»( 1 ) 


-5 


Tew 


V 1 Mpc J 




[mi 



Therefore, 



r(T(rj),L) 



10 34 at the EW transition 
10 14 at the QCD transition . 



(21) 



(22) 



Thus, the amplitude of the large scale magnetic fields turns to be within the range 
necessary to be amplified by dynamo models. 

It must be noticed that eq. (l20l) - (l22ll provide the long wavelength contributions to the 
magnetic field production which dominate for L > 0. lpc. 

The results presented here (and in ref.L19]) arise from quantum loop effects taking into 
account quantum processes that only occur out of equilibrium. These quantum processes 
are classically forbidden. That is, they are virtual processes that take place off-shell. 

In summary, the large scale primordial magnetic fields generated during the scaling 
stage after a phase transition are a plausible mechanism to generate primordial magnetic 
fields which are further amplified by the collapse of protogalaxies and by astrophysical 
dynamos. 
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